Phase transitions as topology changes in configuration space: an exact result 
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£NJ . The phase transition in the mean-field XY model is shown analytically to be related to a 

» I ' topological change in its configuration space. Such a topology change is completely described by 

f> |_ means of Morse theory allowing a computation of the Euler characteristic — of suitable submanifolds 

of configuration space — which shows a sharp discontinuity at the phase transition point, also at finite 
N. The present analytic result provides, with previous work, a new key to a possible connection of 
topological changes in configuration space as the origin of phase transitions in a variety of systems. 

PACS number(s): 05.70.Fh; 02.40.-k; 75.10.Hk 

Phase transitions (PTs) are one the most striking phenomena in nature. They involve sudden qualitative physical 
changes, accompanied by sudden changes in the thermodynamic quantities measured in experiments. From a math- 
ematical point of view, both qualitative and quantitative changes at PTs are conventionally described by the loss of 
! 1 * ' analyticity of the probability measures and of the thermodynamic functions. According to statistical mechanics, in 
the grandcanonical and canonical ensembles, such a nonanalytic behavior can exist only in the thermodynamic limit, 
i.e., in the rather idealized case of a system with N — > oo degrees of freedom jjj. PTs in real systems would then be the 
"shadow" , at finite but large N, of this idealized behavior. The way nonanalytic behavior can emerge as N — > oo due 
to singularities of the probability measures of the statistical ensembles has been rigorously studied by Yang and Lee in 
the grandcanonical ensemble Q and by Ruelle, Sinai and others || in the canonical ensemble; however, to the best of 
i-^j | our knowledge, no equally rigorous approaches to this problem exist in the microcanonical ensemble [Q. Moreover, the 
necessity of taking the N — > oo limit to speak of PTs seems less satisfactory today, since there is growing experimental 
evidence of PT phenomena in systems with small N (e.g., atomic clusters and nuclei, polymers and proteins, nano 
and mesoscopic systems) . There is also a deeper reason why the conventional approach to PTs is not yet completely 
satisfactory. Consider a classical system described by a Hamiltonian Ti. = K(n) + V(ip) where K(ir) — | X^)Li ^1 i s 
the kinetic energy, V(ip) is the potential energy and ip = and 7r = {7ri}'s are, respectively, the canonical conjugate 
| coordinates and momenta. Although in principle all the information on the statistical properties is contained in the 

■ function V(ip), no general result is available to specify which features of V(ip) are necessary and sufficient to entail the 
£SJ ', existence of a PT. This is the more surprising since in many cases, knowing a priori that a system undergoes a PT, 

■ several relevant properties of the PT can be predicted just in terms of very general features of V(<p) (e.g., by means 
of renormalization-group techniques). In the light of some recent results p|-pd| an alternative new approach approach 
seems actually possible by resorting to topological concepts: PTs would then be related to topology changes (TCs) of 
suitable submanifolds of configuration space, defined by the potential energy V(ip). The connection between topology 
of configuration space and PTs can be heuristically understood as follows. The microcanonical volume £l(E) can be 
written, at total energy E, as a convolution integral Jl2| 
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n(E) = J* dK Vol (§^ 1/2 ) Vol (M E -k) , (1) 



where §^ is an ?i-dimcnsional hypersphere of radius R, K represents the kinetic energy and My = Me-k is the 
subset of the configuration space M contained within the equipotential hypersurface of level V = E — K. As K 
is varied, all the hyperspheres are trivially topologically equivalent. On the contrary, the topology of Me-k can 
change for some value of K, also at finite (and even small) N. The larger N, the better Q(E) is approximated by the 
product of volumes Vol (§^^1/2) Vol (M(y)), where (K) and (V) are, respectively, the kinetic and potential energy 

averaged over the microcanonical equilibrium distribution and (K) = E — (V). An abrupt TC of the Me-k can 
yield singular derivatives of £l(E) jl(|[l3]]; if it is persistent with N, such a TC will result in a loss of analyticity of 
the thermodynamic observables, only in the A^ — > 00 limit. However, in the present approach the loss of analyticity 
is due to a deeper primitive topological cause of a PT, which is already present at finite A^ jl4| . 

Although it has been proven that TCs in configuration space are a necessary condition for PTs to occur p0| , it 
is known that not all TCs are related to PTs |TT| , and no general argument is yet at hand to define the sufficient 
conditions under which a TC is actually related to a PT. The study of particular models is then crucial to obtain 
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clues towards the general solution of this problem. The main purpose of this Letter is to give a complete analytical 
characterization of the TCs in the configuration space of the so-called mean-field XY model, which may indicate the 
difference between the TC related to the PT and other TCs. Before we discuss this model, we first summarize a few 
needed facts about topology. 

The TCs we arc referring to are those transformations which map a manifold onto a new one which is not diffco- 
morphic to the previous one, i.e., which cannot be mapped back to it by means of a diffcrcntiable transformation. A 
TC is therefore any transformation which "breaks the fabric" of a manifold: making a hole — without boundary — in a 
sphere tranforms it into a torus, and there is no smooth way to transform a torus back to a sphere. Morse theory pl| 
provides a way of classifying TCs of manifolds, and links global topological properties with local analytical properties 
of smooth functions defined on them, so that they can be used as a practical tool to study their topology. Given a 
(compact) TV-dimensional manifold M and a smooth function / : M t— > M, a point x £ M is called a critical point of 
/ if df = 0, while the value f(x) is called a critical value. A level set / _1 (a) = {x 6 M : f(x) — a} of / is called 
a critical level if a is a critical value of /, i.e., if there is at least one critical point x G f~ 1 (a). The function / is 
called a Morse function on M if its critical points are all nondegenerate, i.e., if the Hessian of / at x has only nonzero 
eigenvalues, so that the critical points x are isolated. We now consider the configuration space of a classical system as 
our manifold M, and the potential energy per particle V(ip) = V(<p)/N as our Morse function. Then the submanifolds 
M v of M whose topology we want to investigate are 

M v = V _1 (-oo,«] = {ip £ M : V(ip) < v} , (2) 

i.e., the same as the My = Me~k defined above (where in M v , V has been rescaled by i.e., v — V/N, in order 
to make the comparison of systems with different N easier). All the submanifolds M v of M, with increasing v, have 
the same topology until a critical level V _1 (v) is crossed. Here the topology of M v changes in a way completely 
determined by the local properties of the Morse function: at any critical level a fc-handle is attached Jig ], where 
k is the index of the critical point, i.e., the number of negative eigenvalues of the Hessian matrix of V. Notice that if 
there are m > 1 critical points on the same critical level, with indices hi, . . . , km, then the TC is made by attaching 
m disjoint handles H^ kl \ . . . ,H^ km '. This way, by increasing v, the full configuration space M can be constructed 
sequentially from the M v . Knowing the index of all the critical points below a given level v, we can obtain exactly 
the Euler characteristic of the manifolds M v , given by 

N 

X (M V ) = £(-1) k n k (M v ) , (3) 

fc=0 

where the Morse number ^ is the number of critical points of V which have index k p8| . The Euler characteristic x 
is a topological invariant: any change in x(M v ) implies a TC in the M v . 

Thus, in order to detect and characterize topological changes in M v we have to find the critical points and the 
critical values of V, which means solving the equations 

^ = 0, i = l,...,N, (4) 

and to compute the indices of all the critical points of V, i.e., the number of negative eigenvalues of its Hessian 

11^ = —^— i.j 1 V. (5) 

In the case of the mean-field XY model, which describes a system of N equally coupled planar classical rotators Jl7| , 
due to the mean-field character of the interactions, such a calculation can be done in a completely analytical way. 
This allows then a new discussion of the relationship between TCs and the PT of this model, whose potential energy 
is 

j N N 

where € [0,27r] is the rotation angle of the i-th rotator and h is an external field. The model describes also a 
planar (XY) Heisenberg system with interactions of equal strength among all the classical spins Sj = (cos (ft, sin tpi). 
We consider only the ferromagnetic case J > 0; for the sake of simplicity, we set J = 1. In the limit h — > 0, the 
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system has a continuous PT, with classical critical exponents, at T c = 1/2, or e c — 3/4, where e = E/N is the total 
energy per particle |l7| . Defining the magnetization vector per particle m = (m x , m y ), where m x = i? Y^iLi C0Sl fij 
m y = ~k ^2iLi sm( / 5 ij the potential energy V can be written as a function of m as: 



V(<p) = V(m x ,m v ) = y(l-m|-m^) - hN m x . (7) 



The range of values of the potential energy per particle, V = jj, is then —h < V < 5 + 

The configuration space M of the model is an A-dimensional torus, being parametrized by the N angles {^i} = 
ipi, . . . ,(p N . We now study the topology of the family of submanifolds M v for this model. First, since TCs of M v can 
occur only at critical points of V, there are no TCs when v > i + 4-, i.e., all the M v 's with v > \ + ^- must be 
diffeomorphic to the whole M, that is, they must be TV-tori. Then one has to find all the solutions of Eqs. ([|), which 
can be rewritten in the form 0,0] 



(m x + h) sin tfi — m y cos (pi = , i = 1, .. . ,N . (8) 

As long as {m x + h) ^ and m y 7^ (m x and m y are both zero only on the level v = 1/2 + h 2 /2), the solutions of 
Eqs. (|^) are all those configurations for which the angles ifi are either or n. These configurations correspond to a 
value of v which depends only on the number of angles which are equal to n, and using Eq. (Q) one obtains 



2n„) , (9) 



where < < N. We have thus shown that as v changes from its minimum —h (corresponding to — 0) to \ 
(corresponding to ii T = -j) the manifolds M v undergo a sequence of topology changes at the N critical values w(n,r) 
given by Eq. (||). There might be a TC also at the last (maximum) critical value v c = | + However, the above 
argument does not prove it, since the critical points of V corresponding to this critical level may be degenerate p9[ , 
so that on this level V would not be a proper Morse function. Then a critical value v c is still a necessary condition 
for the existence of a TC, but it is no longer sufficient. However, as argued in Refs. 0,0, it is just this TC at v c 
which should be related to the thermodynamic PT of the mean-field XY model. For, the temperature T, the energy 
per particle e and the average potential energy per particle u — (V) obey, in the thermodynamic limit, the equation 
2e = T + 2u(T), where we have set Boltzmann's constant equal to 1. Substituting in this equation the values of 
the critical energy per particle and of the critical temperature we get u c = u(T c ) = 1/2; as h — > 0, v c — > ~j so that 
v c = u c . Thus a TC in M occurring at this v c , where v c is independent of N, is connected with the PT in the limit 
TV — > 00, and h — > 0, when indeed thermodynamic PTs are usually defined. 

Let us now prove that a TC at v c actually exists and try to understand why it is different from the other TCs, i.e., 
those occurring at < v < v c . To do that we characterize all the TCs occurring at the critical values < v < v c 
using Morse theory, computing the indices of the critical points of V. At these points, where the angles are either 
or 7r, we can write the Hessian matrix (|j) in the form 



where D is diagonal, D = diag(^), with 



NH = D + B (10) 



Si = (m x + h) cos (pi , (11) 



and the elements of B, bij, can be written in terms of a vector a whose N elements are either 1 or — 1: 

kj = -jjViVj , (12) 

where <7i — +1 (resp. —1) if tpi = (resp. 7r). The matrix B has only one nonzero eigenvalue. This implies that the 
number of negative eigenvalues of H equals the number of negative eigenvalues of I? ±1 |2C[| , so that as N gets large 
we can conveniently approximate the index of the critical point with the number of negative <5's at x. At a given 
critical point, with given n w , the eigenvalues of D are 

Si = m x + h i = 1, . . . , iV — n.^ ; (13a) 

Si = -{m x + h) i = N -n„ + l,...,N , (13b) 
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where the ir-component of the magnetization vector is m x = 1 — so that m x > (resp. < 0) if Ji r < =j (resp. > ^). 
Then, if the external field h is sufficiently small, and denoting by index^?,-) the index of a critical point with given 
we can write 



index(n w ) 
index(nTr) 



TV 



if n n < f 
if n,, > £ 



(14a) 
(14b) 



The number C(n,r) of critical points having a given n^, which is the number of distinct strings of 0's and 7r's of 
length TV having occurrences of 7T, is given by the binomial coefficient C(n^) — ( )■ Thus, at any critical level 
—h < u(n7r) _ 5j where t^n^) is given by Eq. (^), a topological change in M v occurs, which is made up of attaching 
C(n w ) fc-handles, where kfa^) = index(n 7r ) given in Eq. (|l4|). Here as a function of v can be obtained by solving 
Eq. (§), yielding 



int 



l + /i± J/i 2 -2(t> 



(15) 



where int {a} stands for the integer part of a. Eqs. (|T3) and dTl) allow us to write the Morse numbers fik of the 
manifolds M v , for —h<v<k + ^,Ba 

Mfc(«) 



l-9(fc-n^ ) («)) + ©(^ 
TV X 



,(+), 



«)) 



A; = 0,1,..., TV, 



(16) 



where O(z) is the Heaviside theta function. We note that since < ni 

] < f and f + 1 < ni +) < TV, Eq. © 

implies 







Vfc> f , 



i.e., no critical points with index larger than ^ exist as long as v < \ + \- On the other hand, for v > \ 
must be an TV-torus , and for any Morse function on such a manifold one has pi] ] 



k = 0,l,...,N . 



Thus, as \ < v < 
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the manifold is only "half" an iV-torus, and since we know that for v > 



1 -l hi 

2 ' 2 



(17) 

(18) 

is 

/c-handles 



an iV-torus, we conclude that atw = w c = | + -^-a TC must occur, which involves the attaching of 
for each k ranging from -y + 1 to N. This is surely a "big" TC: all of a sudden, "half" an iV-torus becomes a full 
TV-torus. Now we can use Eqs. (^), ( |l5| ) and ( |l6| ) to compute the numerical values of the Euler characteristic of the 
manifolds M v as a function of v: it turns out that x jumps from positive to negative values, so that it is easier to look 
at In Figure [j], log(\x\(M v )) / N is plotted as a function of v for various values of N ranging from 50 to 800. The 
"big" TC occurring at the maximum value v c of V, which corresponds in the thermodynamic limit to the PT, implies 
a discontinuity of \\\, jumping from a big value (0(e ) in our case) to zero, which is the value of x f° r an TV-torus. 

The analytical results we have presented provide a possible hint about what could be the sufficient conditions 
for a TC in configuration space to yield a PT. In the model studied here, the TV TCs which are not related to the 
PT involve the simultaneous attachment of handles which are all of the same type, while that occurring at v c is 
the simultaneous attaching of handles of O(N) different types. Hence we might conjecture that this is a sufficient 
condition for a TC to be in one-to-one correspondence with a thermodynamical PT, also in other models. Further 
confirmation or refinements of this conjecture may provide the basis for a full theory of the origin of PTs based on 
the topology of configuration space as encoded in the potential energy. It is worth noticing that this might be useful 
not only to study PTs in the microcanonical ensemble and/or in finite systems, but also in systems which are difficult 
to study with conventional approaches, like disordered ones or glasses. For glass-forming liquids topological concepts 
have been recently invoked |2^] . 
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FIG. 1. Plot of log(\x\(M v ))/N as a function of v for h = 0.01 and increasing TV = 50,100,200,400,800 (from bottom to top). 
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